Low-speed impact craters in loose granular media 
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We report on craters formed by balls dropped into dry, non-cohesive, granular media. By explicit 
variation of ball density pt, diameter D^, and drop height H, the crater diameter is confirmed to 



scale as the 1/4 power of the energy of the ball at impact: D c ~ (p;,!);, 3 H) 1 ^ 4 . Against expectation, 



a different scaling law is discovered for the crater depth: d ~ (p& 3 ^ 2 Db 2 H) 1 ^ 3 . The scaling with 
properties of the medium is also established. The crater depth has significance for granular mechanics 
in that it relates to the stopping force on the ball. 

PACS numbers: 45.70.-n, 45.70.Cc, 83.80.Fg 
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Sand is fragile: If you set down a ball, no matter how 
gingerly or how roughly, the sand can barely support the 
ball's weight; one slight tap and the ball digs in deeper. 
At the same time, sand can be very strong: If you drop 
the ball from some height, the sand can stop it quickly, 
forming a shallow crater in the process; a higher drop 
height will lead to only a slightly deeper crater. Simi- 
larly intriguing combinations of toughness and fragility 
have led to tremendous recent research activity into the 
physics of granular media in general 0,0- At the scale 
of grain-grain interactions, there are only a few possi- 
ble forces. There can be normal forces perpendicular to 
the contact plane, there can be static and sliding friction 
parallel to the contact plane, and there can be inelastic 
collisions. Since granular packings are random, the nor- 
mal forces are random too but can be correlated over long 
distances in so-called "force chains." Dissipation can be 
either through inelastic collisions or sliding friction. But 
which of these is responsible for taking up the energy of 
an impacting ball? 

Geophysicists have long been interested in the craters 
formed by impact or explosion 0, 0, 13 ■ Since meteorites 
generally strike at non-normal angles but nevertheless 
produce circular craters, it is believed that impacts can 
be likened to explosions. Extensive data for buried ex- 
plosives and high-speed (km/s) ballistics impact indicate 
that the crater diameter D c often scales as a power of en- 
ergy. For example, the exponent is 1/3 when the energy 
is dissipated by plastic flow of the medium throughout a 
volume ~ D c 3 ; it is 1/4 when the energy goes into lifting 
a volume ~ D c 3 by a distance ~ D c against the force of 
gravity. This "gravity-limited" regime was recently ob- 
served in low-speed laboratory impact experiments 
There, steel balls of various diameters Dt were dropped 
into sand from various heights H; the resulting crater 
diameters scaled as D c ~ (-D& 3 H) 1 / 4 . For high-speed 
impacts in loose sand-however, ballistics data support 
D.~ ft 1 /'Di 5 /«fli/«0l 

Here we focus on the depth of craters formed by projec- 
tile impact. By contrast, the diameter has traditionally 
been the primary observable. Prior literature assumes 
that depth is proportional to diameter, as in explaining 



the scaling with energy, but this has not been checked in 
laboratory experiments. There are few actual observa- 
tions, but "simple craters" formed by meteorite impact 
are thought to be parabolic with a depth equal to 1/5 of 
the diameter 0] ( "complex craters" are larger and exhibit 
flat floors or central peaks). The crater depth is also a 
key physical quantity because it relates to the force ex- 
erted by the medium onto the projectile. Specifically, if 
a ball of mass m is dropped from rest and forms a crater 
of depth d, then the average stopping force satisfies 



(F)d = mgH. 



(1) 



Note that H is the total drop distance, equal to the sum 
of initial height above the medium plus the depth of the 
crater. This simple experiment therefore gives a direct 
signature of the dissipation mechanics. We find that the 
depth d can be much less than the drop height H; there- 
fore, to a suprising extent, the dissipation force can far 
exceed the ball weight. We also find that crater depth 
and diameter do not exhibit the same scaling with ball 
density, ball diameter, and total drop distance. This con- 
tradicts common assumption |3J, |4j, |5| . Our observations 
could lead to new insight both into fundamental granu- 
lar forces, and perhaps as well as into geophysical impact 
and explosion cratering. Analogy may also be found for 
impact on metals || U 

Our experimental procedures are as follows. First, ei- 
ther a 1 L or 2 L beaker is filled half way with a gran- 
ular medium (results are independent of beaker size). 
The container is then horizontally swirled and lightly 
tapped in order to level the surface without noticeable 
compaction. Next a ball is fixed in the jaws of a wrench 
mounted to a ringstand. The ball is released with zero 
translational and rotational speed directly above the cen- 
ter of the beaker. It splashes into the granular medium 
and comes to rest with some of its lower portion buried 
and its top fully exposed. This produces a crater with 
a circular rim that extends above the original horizontal 
surface of the medium. As is standard, we define the 
crater diameter by the location of maximum rim height. 
But now, we also measure the crater depth. Useful defini- 
tions might include the depth above which the ball is ex- 
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posed, or the depth below which grains are not disturbed. 
Rather, to make contact with Eq. (JTJ and the average 
stopping force on the ball, we define the crater depth as 
the distance of the bottom of the ball below the initial 
horizontal surface of the medium. This is determined to 
±0.5 mm by measuring the location of the top exposed 
surface of the ball. Before proceeding, two bounds on 
the measured quantities should be noted. First, it is not 
possible to measure the crater diameter if it is less than 
the ball diameter. But even if the diameter is less than 
this minimum, the depth is still easily measured. Second, 
a ball dropped from infinitesimally above the surface of 
the medium can still penetrate to a non-zero depth, more 
so for denser balls. This defines a minimum crater depth, 
and corresponds to a minimum in total drop distance H. 
Our dynamic range for H is 1 — 1000 mm at best, but is 
more typically a factor of one hundred. 

In Fig. ^ we show crater size vs loss in ball energy, 
E = mgH, for a variety of balls dropped into dry, 
monodisperse, 0.2 mm diameter glass beads. When so 
plotted, the diameter data collapse onto a power-law 
curve, D c = (0.24 cm/erg 1 / 4 )^ 1 / 4 . The crater diame- 
ter scaling is thus D c ~ (pbD^H) 1 / 4 . This agrees with 
the gravity-limited scaling observation of Ref. @. Our 
data also extend that work, both in dynamic range and 
in explicit variation of the ball density. However, the 
most striking feature of Fig. ^ is that the crater depth 
data do not show a similar collapse! Apparently, the 
crater depth does not scale with the ball energy. Fur- 
thermore, it appears to follow a different power-law with 
total drop distance, d ~ iJ 1 / 3 . The weighted average of 
the exponents for power-law fits vs H are 0.231 ±0.005 for 
diameter and 0.318 ±0.005 for depth; for individual data 
sets, the average uncertainty is 0.03. Against natural ex- 
pectation and contrary to prior assumption 0, 0, 0, @ , 
the crater depth and diameter are separate length scales 
set by different physics. 

Before further analyzing the data in Fig. ^ two spe- 
cific comparisons should be noted. First, the smaller 
nylon ball and the silicon rubber ball have nearly the 
same density and diameter, but have very different sur- 
face properties: nylon is slick and silicon rubber is tacky. 
But as seen in Fig. ^ the crater depth data for these 
two balls (0 and +, respectively) are indistinguishable. 
Therefore, friction between the ball and grains cannot be 
the dissipation mechanism responsible for stopping the 
ball. Second, the "live" and "dead" balls have nearly the 
same density and diameter, but have very different resti- 
tution coefficients. But as seen in Fig. ^ the depth data 
for these two balls (• and ■. repectively) are indistin- 
guishable. Therefore, no significant energy is transfered 
to internal degrees of freedom of the ball. 

We now explore the full form of the scaling laws 
for crater size. If the ball diameter and drop height 
are the only relevant length scales, then the sim- 
plest dimensionally-correct laws consistent with the H- 
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FIG. 1: The diameter and depth of craters formed by balls 
dropped into 0.2 mm diameter glass beads, as a function of 
total energy loss, E = mgH . Each data point represents the 
experimental result for a single drop height H. Plotting vs 
mgH collapses the diameter data, as seen earlier in Ref. [f|, 
onto a 1/4 power law. By contrast, the depth data does not 
collapse, but follows a 1/3 power law. Symbol, ball type, ball 
density, and ball diameter are as follows: o for hollow plas- 
tic, 0.26 g/cm 3 , 2.54 cm; □ for wood, 0.83 g/cm 3 , 1.59 cm; 
for nylon, 1.10 g/cm 3 , 1.59 cm; x for nylon, 1.10 g/cm 3 , 
2.54 cm; + for silicon rubber, 1.10 g/cm 3 , 1.52 cm; A for 
acrylic, 1.20 g/cm 3 , 1.59 cm; • for live ball, 1.20 g/cm 3 , 
3.82 cm; ■ for dead ball, 1.30 g/cm 3 , 3.82 cm; ♦ for del- 
rin, 1.40 g/cm 3 , 1.59 cm; A for teflon, 2.20 g/cm 3 , 1.59 cm; 
T for ceramic, 3.90 g/cm 3 , 1.27 cm. NB: Data for three 
even denser balls are omitted: stainless steel, 7.90 g/cm 3 , 
2.54 cm; lead, 11.3 g/cm 3 , 1.13 cm; and tungsten carbide, 
16.4 g/cm 3 , 1.91 cm. These produce craters that also scale as 
D c ~ if 1 / 4 and d ~ H 1 ^ 3 ; however, they do not exhibit the 
same D c ~ _B 1//4 collapse as the less dense balls. 



dependence of Fig. ^ would be D c oc D^^H 1 / 4 and 
d cx Z>b 2//3 7J 1 / 3 . To test this, we deduce the dimen- 
sionless proportionality constants for all the size vs H 
data. The results are plotted in Fig. E(a) as a func- 
tion of ball density. For the crater diameter, we find 
D c ~ Pb 1 ^ 4 as expected from the above E 1 / 4 collapse. 
However, this trend is violated by the three densest balls 
(stainless steel, lead, and tungsten carbide), which were 
omitted from Fig.^ For these three, the crater diameter 
is D c ps l.bD^^H 1 / 4 with no apparent dependence on 
ball density. For the crater depth, we have no a priori ex- 
pectation other than that denser balls should penetrate 
deeper. The data in Fig.Efa) are consistent with a simple 
power-law form, d ~ Pb 1 ^ 2 ', the uncertainty in the expo- 
nent is conservatively ±0.05. Note that this holds for 
all balls, even the three densest where D c ~ Pb 1 ^ 4 fails. 
As a final consistency check, the explicit dependence of 
crater size on ball diameter is shown in Fig.^b). Though 
the dynamic range is only a factor of four, the data are 
consistent with expectation, D c ~ D^l 4 and d ~ D 2 ^. 

Next we vary the properties of the granular medium, 
studying craters formed by a 1-inch diameter nylon ball 
dropped into the media listed in Table [U Besides grain 
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FIG. 2: Scaling of crater size with ball density and diameter. 
Each point represents the result of a D c oc H 1 ' 4 or d oc ff 1 ' 3 
power-law fit to the size vs H data sets of Fig. divided by 
the labeled combination of ball diameter or density and drop 
height; error bars are comparable to symbol size. All data are 
for 0.2 mm diameter glass beads. Note that the three densest 
balls do not obey the D c ~ pfc 1 ^ 4 scaling; data for these balls 
are boxed in both (a) and (b). 



size, the other specifications are the (bulk) grain density 
p g , and the angle of repose 8 r . Since identical crater- 
ing is found in 0.2 and 1 mm diameter glass beads, the 
grain size is not an important length scale (as suspected 
already based on the ball diameter and drop height scal- 
ing). Unless the ambient air plays a crucial role, the de- 
pendence on p g must be the reciprocal of the dependence 
on pb. We have no such guess for behavior vs 9 r . But 
since tan 8 r is roughly the coefficient of friction between 
grains, it must play a role. To test all this we divide 
out the expected p g dependence and plot the results vs 
tan 9 r in Fig. [21 The dynamic range is only a factor of 
two for both p g and p = tan6? r , so power-law fits give ex- 
ponents to within only ±0.3. Nevertheless, the data are 
consistent with D ~ (p g /i 2 ) -1 / 4 for the crater diameter 
and d ~ (pgp 2 )^ 1 ^ 2 for the crater depth. Comfortingly, 
the same combination of material properties appears in 
both results, and also leads to a stopping force that is 
proportional to p (next). 

The final crater diameter and depth laws established 
by Figs. 11131 are, respectively: 

D c = Qm^/ipgp 2 )] 1 '^^'^ 1 ^, (2) 

d = 0T6[p 6 /( PsM 2 )]^A, 2/3 tf 1/3 . (3) 



Whereas the depth law holds for all our observations, the 
diameter law fails for dense balls, pb > 4 g/cm 3 (at least 



TABLE I: Specifications of the granular materials. The den- 
sity refers to the bulk material, not the individual grains. The 
angle of repose, r , was measured by the draining method |Tcj|. 
which seems most appropriate for craters; it gives the grain- 
grain friction coefficient as p — tanf? r . The beach sand is the 
only significantly polydisperse material. 
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FIG. 3: Scaling of crater size with grain-grain friction coeffi- 
cient, p = tanf? r , for l-inch diameter nylon ball dropped into 
the media specified in Table |l] Each point represents the re- 
sult of a D c oc H 1 / 4 or d oc H 1 ^ 3 power-law fit to the size vs 
H data sets, divided by the labeled combination of ball and 
grain densities, ball diameter, and drop height; error bars are 
comparable to symbol size. 



for glass beads - the materials dependence of this break- 
down has not been investigated). For a fixed granular 
material, the diameter data collapse onto a 1/4 power- 
law when plotted vs pbDj^H (as in Fig.f^); therefore, the 
diameter scales with ball energy at impact. By contrast, 
the depth data collapse onto a 1/3 power-law when plot- 
ted vs pb 3 / 2 Db 2 H (not shown - but since the scatter for 
depth and diameter is the same in Figs. I2l.il the qual- 
ity of collapse is also comparable); therefore, the depth 
scales as neither ball energy nor ball momentum at im- 
pact. The diameter and depth are separate lengths set 
by separate physics. 

The crater depth d may be the more fundamental 
length scale, in that it naturally relates to the under- 
lying granular mechanics. Energy conservation, Eq. 
and the depth law, Eq. ©, give the average stopping 
force acting on the ball as 



mg 



6.3„ ( *L 
\Pb 



1/2 / \ 2/3 



(4) 



Evidently, it can be much greater than the ball weight, 
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mg, which the granular medium can barely support in a 
static situation. These results help constrain the form of 
the instantaneous stopping force on the ball throughout 
the impact process. If it varies only with ball speed, 
independent of depth, then it must scale as F ~ v 4 / 3 . 
Alternatively, if the stopping force varies only with ball 
depth, independent of speed, then it must scale as F ~ 
z 2 . Solution of F = ma gives d ~ H 1 / 3 for both cases. 

We now discuss these possible force laws in the con- 
text of the current understanding of granular mechanics. 
First, the kinetic theory of granular h ydr odynamics gives 
a rate-dependent stopping force [Tl], [l^] ■ In Ref. it 
was used to analyze the peak height of the granular jet 
formed when a heavy ball strikes the medium and be- 
comes deeply submerged. There, the viscosity of the 
medium is proportional to the shear rate, which scales 
as ball speed divided by diameter. This gives a stopping 
force of F a p g D g 2 v 2 , where D g is the grain diame- 
ter and the numerical constant is set by the grain-grain 
restitution coefficient. This cannot account for our ob- 
servations, since the dependencies on grain density, grain 
size, ball size, and ball speed are all wrong. Furthermore, 
a F ~ v 2 drag force is not even strong enough to bring 
an object to rest! A modification, whereby the viscosity 
increases with packing fraction |14j . would be required. 

Second, the lateral drag force on an object slowly 
pulled through a granular medium at constant depth 
is rate- independent, as in plowing a field. Recently it 
was found to scale as the product of the object's cross 
section and the hydrostatic pressure at that depth [lfij . 
If this applies to our work, where the ball moves down 
rather than sideways, then the instantaneous drag force 
is proportional to the weight of displaced grains: F = 
np g giT(DbZ 2 /2 — z 3 /3). Here z is the depth of the bot- 
tom of the ball, the volume of displaced grains equals the 
submerged volume of the ball, and n is a materials param- 
eter like p. For z <C -D&, the leading term is F ~ z 2 and 
the final crater depth is d = [{p b / p g )D 2 H /rj\ 1/3 . This 
is quite similar to our observation, Eq. (J3J; however the 
ball- and grain density dependencies are incorrect. The 
parameter range in Fig. Efa) is great enough to easily 
rule out d ~ pb 1 ^ 3 in favor of d ~ pb 1 ^ 2 - 

Altogether, it appears that our observation for the scal- 
ing of crater depth, Eq. (JSJ, cannot be explained us- 
ing prior work. Friction at the ball surface, collisional 
granular hydrodynamics, and the plowing of hydrostatic 
grains, are all ruled out. One possible scenario is that 
as the ball crashes into the medium, it jams together 
the grains underneath. The normal force between these 
grains thus becomes much greater than the hydrostatic 
pressure. As the ball moves, the grain contacts slide so 
that each dissipates a total amount of energy given by the 
normal force times grain size. New contacts are formed 



as the old ones break. This loading and breaking of force 
chains gives rise to the dissipation force that ultimately 
stops the ball. Another possibility is that dissipation is 
due to sliding friction between force chains and the sur- 
rounding unloaded grains. Further work is needed to 
model either effect. Perhaps the scaling of crater diame- 
ter, Eq. could then be explained as a consequence of 
the stopping force between ball and grains. 

We thank R.P. Behringer, S.R. Nagel, and J. A. Rud- 
nick for helpful suggestions. This material is based upon 
work supported by the National Science Foundation un- 
der Grant No. 0070329. 

Note added in proof. A similar experiment has been 
reported, in which a steel ball was dropped into glass 
beads of several different sizes jj^. There, the impact 
energy was greater and the granular packing may have 
been looser, so that the ball always became submerged. 
The crater diameter scaled as iJ 1 / 4 ; the crater depth 
(measured to the bottom of the crater, not the bottom 
of the ball) was about 1/8 the diameter but increased 
with H. The crater morphology and the analogy with 
planetary craters were emphasized. 
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